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Abstract 

I We study a model of scalar field with kinetic terms non- minimally coupled to to 

I the curvature, and the scalar field coupled to the Gauss Bonnet 4-dimensional 

invariant. A variety of solutions are found, giving rise to different cosmologi- 
cal scenarios. A unified description of early time matter (radiation) dominance 
I with transitions to late time quintessence and phantom phases have been made. 

Among others, we found solutions unifying asymptotically the early power law 
behavior and late time cosmological constant. Solutions of Chaplygin gas and 
IL" generalized Chaplygin gas cosmologies have also been found. 
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1 Introduction 

The late time acceleration of the universe is one of the most important problems of 
modern cosmology, which is supported by astrophysical data from distant la super- 
novae observations [1], [2], [3], cosmic microwave background anisotropy [5j, and 
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large scale galaxy surveys [6]. The interpretation of astrophysical observations indi- 
cates that this accelerated expansion is due to some kind of negative-pressure form 
of matter known as dark energy ( [7] , [8] , [9] ) . The combined analysis of cosmological 
observations also suggests that the universe is spatially flat, and consists of about 
~ 1/3 of dark matter, and ~ 2/3 of homogeneously distributed dark energy with neg- 
ative pressure. The dark energy may consist of cosmological constant, conventionally 
associated with the energy of the vacuum [10], [H], or alternatively, could came from 
a dynamical varying scalar field at late times which also account for the missing en- 
ergy density in the universe. A widely explored scalar field models are quintessence 
fr2\ . [13], [S], tachyon [151, phantom [T^, K-essence [IT], [IB], and dilaton [TP] (for a 
review see [7j). The scalar fields are allowed from several theories in particle physics 
and in multidimensional gravity, like Kaluza-Klein theory. String Theory or Super- 
gravity, in which the scalar field appears in a natural way. Another alternative to the 
explanation of the DE is represented by the scalar-tensor theories, which contain a 
direct coupling of the scalar field to the curvature, providing in principle a mechanism 
to evade the coincidence problem, and naturally allowing (in some cases) the crossing 
of the phantom barrier [20], [21]. From a pure geometrical point of view, the modified 
gravity theories, which are generalizations of the general relativity, have been widely 
considered to describe the early-time inflation and late-time acceleration, without the 
introduction of any other dark component, and represent an important alternative to 
explain the dark energy (for review see [22] and references therein). 
In the present work, we consider a model with non-minimal coupling to gravity 
[23] EH |25], and with additional Gauss Bonnet coupling, specifically we focus in 
a scalar field model with kinetic term non-minimally coupled to gravity and to itself 
[2S1 [23 12H] , with a new term containing the Gauss Bonnet (GB) 4-dimensional invari- 
ant coupled to the scalar field, with the coupling given by an arbitrary function of the 
field. Despite the fact that the GB term is topologically invariant in four dimensions, 
and hence, by itself does not contribute to the equations of motion, nevertheless it 
affects cosmological dynamics when it is coupled to a dynamically evolving scalar 
field. Besides this, if the GB term is coupled to the scalar field through arbitrary 
function /(0), then this is the unique quadratic combination of the Riemann curva- 
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ture tensor, that does not increase the differential order of the equations of motion 
(i.e. GB produces only terms which are second derivatives of the metric in the field 
equations). Therefore, the coupled GB term which preserves the theory ghost free, 
seems as a natural generalization of the scalar field with non-minimal kinetic coupling 
to curvature, which under certain condition is a second order scalar tensor theory. 
Both couplings my have origin in the low energy limit of higher dimensional theories. 
Thus, the kinetic coupling appears as part of the Weyl anomaly in = 4 conformal 
supergravity [291 ED] , while the GB coupling arises naturally in the leading order of 
the a' expansion of string theory [31], [32|. Besides that, the kinetic couplings to 
curvature are also present as quantum corrections to Brans-Dicke theory [33] and in 
different frames in modified gravity [22j . 

Some late time cosmological aspects of scalar field model with derivative couplings to 
curvature have been considered in [31], [35], [36]. On the other hand, the GB invariant 
coupled to scalar field have been extensively studied. In [37] the GB correction was 
proposed to study the dynamics of dark energy, where it was found that quintessence 
or phantom phase may occur in the late time universe. Accelerating cosmologies with 
GB correction in four and higher dimensions have been discussed in [38], [39], [iO] . 
The modified GB theory applied to dark energy have been suggested in [H], and 
different aspects of the modified GB model applied to late time acceleration, have 
been considered among others, in p2], [12], [13], [H]. 

All these studies demonstrate that it is quite plausible that some scalar-tensor cou- 
plings predicted by the fundamental theory may become important at current, low- 
curvature universe. In the present paper we study different cosmological solutions 
of the model with Gauss Bonnet coupling to the scalar field and kinetic coupling 
to curvature. Different solutions unifying early time decelerated behavior, late time 
accelerated quintessence and phantom phases will be studied. In particular, solutions 
describing the Chaplygin gas and the generalized Chaplygin gas have been considered. 
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2 Field Equations 



Let us start with the action, for scalar field kinetic terms non-minimally coupled to 
curvature and coupled to Gauss Bonnet (GB) curvature 



S= [ d'x^\-^R - ld,<pdy - l^R {F,{cl>)d,m) - 

J LlbTTG 2 2 ^2 1) 



where G is the 4-dimensional GB invariant Q = R^ - ARf.^R^"' + R^yp^R^^'P" , S, 



is 



the dark matter action which describes a fiuid with barotropic equation of state. The 
dimensionality of the coupling constants C, and rj depends on the type of function 
Fi{(f)), and the coupling -^2(0) is dimensionless. Besides the couplings of curvatures 
with kinetic terms, one may expect that the presence of GB coupling term may be 
relevant for the explanation of dark energy phenomena. 



Taking the variation of action (2.1) with respect to the metric, we obtain a general 
expression of the form 

R,.-l9,.R = n'[T;:i + T,,] (2.2) 

where = SttG, TJJl is the usual energy-momentum tensor for matter component, 
the tensor T^^^, represents the variation of the terms which depend on the scalar field 
and can be written as 

T — T't' -L -L -X- t'^^ (9 

where T^^, correspond to the variations of the standard minimally coupled terms, 
T^^, TJl^ come from the ^ and f] kinetic couplings respectively, and comes from 
the variation of the coupling with GB. Due to the kinetic coupling with curvature 
and the GB coupling, the quantities derived from this energy-momentum tensors will 
be considered as effective ones. The variations are given by 

= V^0V.0 - lg,uVx<pV'<p - g^.V{<p) (2.4) 



T^. =e [ (^R,. - Ig^uR^ (Fi(0) Va0V^0) + g,.VxV' (Fi(0) V>V^ 
- ^(V^V. + V.V^) (Fi(0)Va0VV) +^(^i(0)V;,</.V.0 



(2.5) 
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- I (VaV^ (Fi(0)VVV.0) + VaV. (Fi(0)VVV^</.)) 
+ ^VaV" (Fi(0)V^0V,0) + ^^7^.VaV^ (Fi(0)V"0V^ 



(2.6) 



and 

t;:^ =4([V^V,F2(0)]i? - (7^.[VpV''F2(0)]i? - 2[VPV,F,{(f))]R,, - 2[V''V,F2(0)]i?,, 
+ 2[VpVPF2{(f))]R^, + 2g^,[V''V''F2{(j))]Rpa - 2[V''V'^F2(0)]/?^p,^) 

(2.7) 

In this last expression the properties of the 4-dimensional GB invariant have been 
used (see |37], [IS]). Variating with respect to the scalar field gives the equation of 
motion 



- -^d^ (ei?Fi(0)5^0 + r^i^'^-Fi (0)9,0 + 9^0)] + ^+ 

H F HP 

— i {^Rdpcpd^cp + vR^^d^-^d"^) - -^g = 

00 00 

Considering the spatially-fiat Friedmann-Robertson- Walker (FRW) metric, 



(2.8) 



ds' = -dt' + a{ty {dr^ + r^dn') 



(2.9) 



And assuming an homogeneous time-depending scalar field , the (00) and (11) 



components of the Eq. (2.2 ), from (2.3 2.7) take the form (with the Hubble parameter 
H = a/a) 



= ^Peff 



with pe// given by 



rl 



Peff 



02 + r(0) + %H^F^{(t))^^ + 3(2^ + r])HFi 



3 dF dF 

- 3(2^ + V)HF,{<P)<P4> - 2(2^ + V)H^^' - 24i/'^0 



(2.10) 



(2.11) 



and 



-2H - 3H^ = K^Peff 



(2.12) 
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with peff given by 



Peff 



12 

+ 4(e + v)HF, 



+ (2e + r/) 



(2.13) 



where we have assumed scalar field dominance (i.e. TJ^ = 0). The equation of motion 



for the scalar field (2.8) takes the form 



• dV ... 

+ 3i/0 + — + 3(2^ + ri)HFi{(P)(l) + 3(14^ + 5r])HHFi 



+ ^(2e + r/)if (2FM4>+^<P' 



1{A^ + V)H' (2Fi(0)0+^02 



(2.14) 



+ 9(4^ + v)H^Fi 



24 HH^ + H 



2 , tt4 \ dF2 







where the first three terms correspond to the minimally coupled field. In what follows 
we study the cosmological consequences of this equations, under some conditions that 
simplify the search for solutions. 



An important simplification of the Eqs. (2.10 2.14) takes place, under the restriction 
on ^ and 7] given by 

r/ + 2^ = (2.15) 



Under this restriction all the kinetic couplings that appear in the action (2.1) become 



summarized in the term G ^^d^ (pd^ (p , where G^y = R^u — \g^yR. In this case the field 



equations (2.10 2.14) contain only second derivatives of the metric and the scalar field, 
avoiding problems with higher order derivatives [M] . The modified Friedmann 



equations (2.10) and (2.12) take the form 



K 



dF2 ■ 



(2.16) 
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and 



-2H - m'^ =K^ 



1 . 



2^H 2Fi 



+ 



o dFo ■■ o (P F' 



The equation of motion reduces to 

■ dV , , 

<p + 3iJ0 + ^ + ?>^H^ 2F 



dF V 



l2^HHF{(p)^ - 24 ( ^ ^ 



+ 
dq 

dF2 



+ 18^H^F{4>),p+ 



(2.17) 



(2.18) 







Note that, independently of Fi{(f)), and for ^2(0) = const., assuming an asymptotic 
behavior of the scalar field as = 0o = const., gives rise to de Sitter solution, as 



can be seen from Eqs. (2.10) and (2.14). From (2.10) and (2.14) it follows that 
V = Vo = const and H = Hq = ka/Vo/3. 



In the following sections we study cosmological solutions of Eqs. (2.16) and (2.18), 
giving rise to accelerated expansion and describing late time cosmological dynamics, 
according to current observations. 



3 Cosmological solutions in the time variable 



We start studying some solutions to Eqs. (2.16) and ( |2.18 ), in the important case 
when the scalar field potential is absent (i.e. V = Q). This leaves the two couplings 
Fi{(f)) and -^2(0) as the degrees of freedom that may characterize the cosmological 
dynamics. Predicting this couplings is very difficult task, but at least for the GB 
coupling there are some indications from fundamental theories like strings, that favor 
the exponential coupling (dilaton) jl6] . Making 0^ 
it follows 

JF2 



1 ih 
~ 3^ I8^m 3^ dt 



i), from Eq. ( |2.16[ ) with V = Q, 

(3.1) 



replacing in (2.18) after simplifications, we obtain 



.dijj 



dH 



,dH dF2 



+ GH^t/j -tlJ^ + 48H^- - 

dt dt dt dt 



+ 72H^ 



dt 



+ 2AH' 



(P F') r,dH J, 

—^ + 12H^— + 18H^ 
dt'' dt 







(3.2) 
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here we have set = 1. Let's consider some solutions to this equation. In the 
following solutions we will consider the GB coupling 



dt H{t) 



(3.3) 



which as will be seen, permits to find some interesting solutions. Replacing (3.3) in 



(3.2 ) one obtains 



h'^ + 6i/V - + 12(2(7 + l)H^'^ + 18(4(7 + l)H^ = (3.4) 

dt dt ^ ' dt ^ ^ > V ; 

Bellow we study some viable late time cosmological solutions that can be extracted 
from this equation. 
1) Constant EoS. 

Here we consider the following behavior for the kinetic term 

i){t) = \^H{tf (3.5) 



replacing in (3.4) one finds the solution 



^ + 24(7 + 12 

^(') = ^^' ^= 6(^ + 12(7 + 3) ^'-'^ 

where C is the integration constant. This solution gives a constant effective equation 
of state (EoS) Weff = —1 + If ?7 > 0,C < 0, then H{t) does not present time 
singularities, presenting power law behavior at late times. Nevertheless the EoS is 
constant all the time, and there is not transition deceleration-acceleration. If (77 < 
0,C < 0), then the solution describes pure phantom phase with Big Rip singularity 



at to = j37], [IE]. The scalar field can be found from (3.5) as 



= Ar/ln(t-Cr7) + 0o (3.7) 
this expression allows to write the time in terms of the scalar field as: t = e^^^^^^ +Cr] 



(making 0o = 0). The GB coupling in terms of the scalar field is obtained using (3.3) 



F2{(f)) = l^e'*/^^^) + - l)e*/(^'') (3.8) 
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where we neglected the constant term which gives tribial (constant) couphng. The 



kinetic couphng is found from (3.1) 



SO, if we take = 6, then Fi becomes constant given by F2 = 8g/{3C,). 
2) Power law solution. 

In this known case the Hubble function is given by 

H{t) = ^ (3.10) 
which gives accelerated expansion provided p > I, and also has constant EoS (wg// = 



— 1 + 2/(3p)). Replacing in (3.4) one finds the solution 



this expression takes a simpler form if we set the integration constant C = and 
g = —3/4: iplt) = 6p^/t^. In this case the scalar field becomes 

(p = VQplnt + (po (3.12) 

This solution allows to express the cosmological time t in terms of the scalar field as 







t = e^p (taking 0o = 0). The GB coupling from (3.3) becomes 



3 ^ 

F2i<P) = -— ev^p (3.13) 
op 

Note that although F2 increases with time {F2 = gt'^/{2p)), the GB term decreases 
faster that F2 increases, giving adequate behavior at late times. The kinetic coupling 



from (3.1) becomes constant, given by Fi = —2/^. 
3) Time dependent EoS I. 

An interesting solution which provides the transition from the decelerating universe 
to the accelerating one, is (see [22] ) 

H{t)=^ + ^ (3.14) 
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with positive constants (7,p). This Hubble function describes early time (t — i- 0) 
power law behavior a cct^ and late time dark energy dominated by the cosmological 
constant H The corresponding effective EoS is 

2p 



w 



-1 + 



(3.15) 



3(7t + p)2 

Here and in what follows, we use the symbol w to mean Weff- The transition from 
decelerated to accelerated expansion occurs at tc = (i/p ~ p)/l- Note that in order 
for this transition to occur, p should be in the interval < p < 1. For a given current 
EoS Wq we can find a relation between and the current time to (age of the universe). 



Using the above expression for tc and (3.15) at to we find 



3v/rT^(i-^) 



(3.16) 



- 3-^1 + Wo\/P 

From this expression we can see that the closer to —1 wq is, the earlier the transition 



occurs. Replacing in (3.4) and solving, one obtains the kinetic term ip. But the 
solution takes simpler form if we choose g = —3/4, and set the integration constant 
equal to zero, which gives 

= ^Jl±Jll (3.17) 



This equation can be easily integrated to find the scalar field (ignoring the integration 

constant) 



(j)[t) = Ve-ft + 



t 



Solving this Eq. with respect to t, one finds: 



t 



p 

7 



W 



0/(V6p) 



(3.18) 



(3.19) 



where W is the Lambert VT-function. The GB coupling, from (3.3) in terms of the 
scalar field is given by 

3p 



4^2 



In 



pil + W 



</>/{V6p) 



w 



</./(V6p) 



(3.20) 



and the kinetic coupling F\ for this solution becomes constant F\ 
4) Time dependent EoS II. 
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Let's consider the following evolutionary scenario 

H{t) 



p _j_ 
t - 



t 



(3.21^ 



The first term guarantees the power law behavior at early times a oct'^, and the second 
term despite the fact that gives rise to Big Rip singularity ait = ts, has the advantage 
that provides the crossing of the phantom barrier (see [22] )• Therefore, this model 
presents three phases: early time decelerated {p < 1), transition from decelerated to 



quintessence phase, and transition from quintessence to phantom phase. From (3.21) 

tp 



the scale factor is a(t) = ag 



w 



= -1 



-, and the EoS is given by 
, 2tsp{ts-2t) + ip~q)t^ 



(3.22) 



3 {tsP +{q- p)tf 

we can obtain a relationship between the different critical times as follows. If we 
consider that at current time to the EoS reaches the value Wq = —1 (which could be 



true according to observational data), then from (3.22) we find (in fact there are two 
roots, but we choose one) 



and the transition decelerated-accelerated expansion occurs at t 



(3.23) 



tc, when Wc 



-1/3. From (3.22) we find 



{p{p -q-l) + ^/pq{q ~ p + 1)) 
{q ~ p){q - p + I) 



(3.24) 



from (3.23) and (3.24) we find the relation between the age of the universe to and t^. 



to 

tr 



(3.25) 



Thus, if we take p = 2/3 which simulates matter dominance at early times, then 
for q = 1/8 one finds to/^c ~ 1-6, which is an appropriate relation according to 
observations. If we start from radiation dominance p = 1/2, then for g 1/4 one 
also finds to/^c ~ 1-6. On the other hand, from (3.23) it follows that the Big Rip 
singularity occurs at the future and depends on the fraction q/p (Big Rip singularities 
in modified F{G) gravity have been studied in |1S]). At this time the EoS reaches 
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the value Wg = —1 — 2/(3g). Hence, the unification of initial decelerated, accelerated 
and phantom universe may emerge in this model. Let's turn to the scalar model that 
gives rise to this scenario. 



The kinetic term ip is obtained by replacing (|3.21|) in (|3.4|). The resulting equation 
has a simpler solution if we take g 
{q - p)t + pts 



m 



-3/4, which is 

[{Qpt^^ + C{ts - tfH) ts - (6(p - q)& + C{ts - tfH) t] 



(3.26) 



One can make one more simplification if taking C = 0, giving 

Considering this last case, we can explicitly integrate to obtain the scalar field as 



76 In 



its - ty 



The GB coupling from (|3.21|) is 
3 



F2{t) = - 



qts 



+ 



pqtl 



4 \{p-qY 2{p-q) {p - qf 



In [p{ts -t) + qt] 



(3.28) 



(3.29) 



and the kinetic coupling from (3.1) and (3.27) is found again as the constant Fi = 
—2/^. In the special case q = p (see |17]) we can express the GB coupling exphcitly 
as function of the scalar field. In this case from (3.28) one can express the time as 
(setting q = p) 

t 



1 + e-«^/(v^p) 



(3.30) 



and integrating the Eq. (3.3) with p = q one finds the GB coupling as 

^2 g20/{V6p)(^g<^/(V6p) _^ 3^ 



i^2(0) 



Bp (1 + e'f'/^^P^y 



(3.31) 



Note that the particular solution (3.27) for the kinetic term is actually ip = 6H . 



4 Cosmological solutions in the e-folding variable 



In this section we study some solutions to Eqs. (2.16) and (2.18) in the x = Ina vari- 



able. It is of interest to derive solutions that not only reproduce the current acceler- 
ated expansion according to observations, but that also have appropriate asymptotic 



12 



behavior at early times. In the e-folding variable x = \na, we have d/dt = Hd/dx. 
Then (p^ = H^{d(j)/dxf and the Eqs. ( |2.16[ ) and ( |2.18| ) take the form 

,2 



and 



1 d 



K 

y 



1 d P 

-H^e + V + %H^Fie - 2AH^^ 

2 dx 



(4.1) 



2 (ix 



12H' 



,dH'^ dF2 

diOC dec 



,dF2 
dx 







where we have multiplied the Eq. (2.18) by and represented {d(j)/dx) 



(4.2) 



e{x). 



Replacing the product Fi9 from (4.1) into (4.2) we arrive at the following equation 

dx g-j 



dx dx 

.dH^ dF2 



AH' 



,dV 



36 

Z2' 



dx 
•AF2 



2V- 



dx 



+ —H^ + 72H^— + UAH^^ + A8H 



:d^F, 







dx dx ' dx ' dx'^ 

This is a first order differential equation for the functions if, 9 and V , and of second 
order for the coupling F2 (in fact is a first order for dF2/dx). To have a solution 
we need to give three of this functions and integrate the equation with respect to 
the remaining one. It is desirable to have some definite criteria to choose the correct 
couplings and the potential in order to prove the predictive power of the theory. There 
are some indications from the fundamental theories like superstring and supergravity 
that suggest potentials and GB coupling as exponentials of the scalar field, but strictly 
speaking, we still lack a definitive criterion for selecting potentials and couplings in 
the scalar tensor theory. Perhaps a good criterion to select the kinetic coupling as 
inverse of the squared scalar field (Fi = 1/0^) [22] is that the coupling constant 
^ becomes dimensionless, but though this is the simplest, is not the unique choice. 
As we have seen in the time variable, a simple exponential or even constant kinetic 
coupling, led to interesting cosmological solutions. Concerning the GB coupling, note 



that in (3.8) and (3.13) the GB coupling is expressed as exponential of the scalar field. 



while in (3.20) and (3.31) is expressed through exponentials of the scalar field. As 
we will see bellow, the choice of the GB coupling, allows to describe the cosmological 
dynamics compatible with current observational data. To begin, we choose the GB 
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coupling satisfying the restriction 

dF2 g 



(4.4) 



dx 

where g is a constant of dimension L^^. This couphng lowers the power of H by two 



and four in the last three terms of Eq. (4.3), making it more tractable. Replacing 



(4.4) in (4.3) one obtains 



dx dx dx dx f a r\ 

Kj dx Kj dx 
Let's now consider different cosmological scenarios that we can extract from this 
equation (after solving this equation we can find the kinetic coupling Fi from Eq. 



(4.1). In this paper we consider an important class of solutions obtained in absence 
of potential. 

In what follows we will introduce the scaled variables, but keeping the same symbols, 
namely we consider the Hubble function as the scaled one, i.e. H — )■ H/Hq, the 
scaled potential energy V — > k'^V/Hq, and the scaled kinetic term as will be clear 
bellow 6 — )■ k'^O. Keeping the same symbols for these functions is also equivalent to 
set = 1 and Hq = 1, where Hq is the current value of the Hubble parameter. After 



this rescalings and making = in (4.5) it follows 



2H^f. + 9^ + ueH^ + 12(1 + 2g)^ + 36(1 + Ag)H^ = (4.6) 
dx dx dx 

This equation depends on two unknown functions, and we have to propose or constrain 
one of them in order to integrate the remaining one. 



Solutions for a given 9{x] 
1) Constant kinetic term. 



Let's consider 6 = const. = X . Replacing in (4.6), we find 



rr2 ^ -2ax 6 12^ + + 3 , ^ ^, 

= Ce cr = ^ — (4.7) 

24^ + A2 + 12 ^ ' 
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This solution corresponds to power law behavior a{t) oc t^^", describing presureless 
matter for cr = 3/2, radiation for a = 2, and accelerated expansion provided 1/cr > 1. 
From the expression for 6, it follows for the scalar field (up to additive constant): 



= Ax. Integrating (4.4), the GB coupling becomes 

m) = §e^- = ^e^5^ (4.8) 
here we neglected the constant of integration, since a constant coupled to GB term 



gives trivial results, due to the geometrical properties of the GB term. From (4.1) we 
find the expression for the coupling Fi 



2a , 



ZCi \ \^ 6 

so, in absence of potential, exponential couplings give rise to power law behavior. 
2) Exponential behavior. 

Let's consider the exponential function for the kinetic term 

e{x) = A^e""^ (4.10) 



replacing in Eq. (4.6) we find the following solution for the Hubble function 

^2 _ f Aj'-^- j_ r,Ja+2i^).V _ 9 - 4g{a -3) -2a 



An interesting case of this solution takes place when 7 = 1, giving 

H{xy = Ae2("-6)^ + Se(3"-^2)^ (4.12) 

taking a = 4 gives a solution that describes two asymptotically important cases: 
radiation dominance at high redshifts (x — )■ — oo) and cosmological constant at far 
future. The value a = 9/2 describes matter dominance at high redshift and phantom 
behavior (negative power law) at far future (if 7 = 1). In this case if we take A = 0.3, 
then the transition decelerated-accelerated expansion occurs at 2;t ~ 0.59, with current 
EoS Wq = —1.05, presenting quintom behavior (see bellow). The equation of state 



that follows from solution (4.11) with the flatness condition A + B = 1, is (in the 
redshift variable, using 1 + z = e~^) 

w(z] = 3 a 777 — — 4.13 

^ ^ 3 3(l + A[(l + z)"-l]) ^ ' 
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the important limits of the EoS are for positive a. The asymptotical EoS at early 
times is 

lim w(^z^ = Wc 

z— >oo 

while 



2 

3--a 



-a 



07 

3" T 

describes the asymptotic future value of the EoS. The current EoS is 



lim w{z) 

2— !>— 1 



(4.14) 
(4.15) 



w{z = 0) = Wq 



2 1. 
:a- - 1 



A)a'j 



(4.16) 



3 3 

For a given set of parameters {a, 7, A), if the solution starts in a decelerating regime, 
the transition deceleration-acceleration occurs at the redshift 



-H-(l-A)i/° 



10 — 2a — «7 



l/a 



(4.17) 



2A{a - 5) 

Thus for instance, we can reach similar results as obtained with the ACDM model: 
if we take a = 9/2, then Woo = 0, describing the matter dominance era, and tak- 
ing 7 = 2/3 the universe evolves towards cosmological constant phase (de Sitter) 
w_i = —1. Taking 7 = 1 gives w^i = —3/2 and the solution evolves from mat- 
ter dominance, passing through quintessence, to finally enter in the phantom phase. 



presenting quintom behavior. Note that in principle, using Eqs. (4. 14)- (4. 16), we 
can satisfy any given in advance asymptotical (1^00,^^-1) and current (wq) behavior, 
within observational expectations. If we remove the GB coupling (i.e. g = 0), then 
the power 7 can not be chosen independently of a, which makes the solution not 
interesting. 



Integrating in (4.10) one finds the scalar field 

2A 



a 



and from (4.4), the GB coupling is 

g 4(0-6) 



2(6 - a)A-^ 



axj2 



12 

-2 + -,7, 
a 



12 

a 



B 



where we have used (4.18), setting 



0, and rescaled 



coupling is found from Eq. (4.1) as 



, 4(a-6) 



Fi{4>) 



+27 



3^ (A02 + Bf \ A20- 



1 + 8(7 



(4.18) 

(4.19) 
The kinetic 

(4.20) 
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Solutions for a given H{x) 

The cosmological equations can also be consistently solved by a given behavior of 
Hubble function (or the scale parameter). In this case we reconstruct the scalar field 
and the couplings. Let's consider some important solutions 



1) ACDM-type solution. 
The Hubble function 



(4.21) 



reproduces the power law behavior at early times, describing radiation or matter 
dominance, and late time behavior governed by the cosmological constant. It can 
also serves as purely dark energy model for appropriate values of a. Replacing in Eq. 



(4.5) and solving with respect to 9 one obtains 



e{x) =C (^]l + fiae"^)"^^^ e(""^^)^/=^ - 3(1 + 4g)- 



3{3 + 4g)aV^i 



-1/2 



116 16 1^26° 



12 -a 

this solution is considerably simplified for g = —3/4, giving 

e{x) =C{Qi + 1^26"^)"^/^ e("-^2)x/2 ^ g 



(4.22) 



(4.23) 



and even simpler if we take the integration constant C = 0. In this case 6{x) = 6, 
and after integration the scalar field is obtained as = ^/6x + (po- In fact is easy to 



check directly, that this particular value 6 = Q satisfies the Eq. (4.5) with H given 



by (4.16). 



For this particular solution, the scalar field is given by = \/6x + (po, and from (4.4) 



and (4.21) we may write the GB coupling explicitly in terms of the scalar field, as 

3 



In 



from Eq. (4.1) we find the kinetic coupling 



f2i + il2e 



(4.24) 



3^ Vli + n2e"<^/^ 



(4.25) 
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Therefore, even without the cosmological constant nor cold dark matter in the action, 
the cosmology of the ACDM model can be reproduced by the combined effect of the 
kinetic and GB coupling. 
2) The Chaplygin gas solution. 

The Hubble function for the Chaplygin gas has the known form [50] (setting = 1 
and with appropriate normalization of A and B) 

= [A + Be-^^Y^^ (4.26) 

This solution comes from the EoS obeyed by the Chaplygin gas, p = And has 
the known properties of describing the presureless matter dominance epoch of the 
universe at early times (a << 1, normalizing the current value of a to 1), and the 
future universe dominated by the cosmological constant, evolving towards de Sitter 
phase at a >> 1. 



Replacing (4.26) in (4.5), we arrive at the equation 

2{A + 5e-^^) — + 3{4A + 3Be-^'') e{x) + 36 (A + 4Ag + 25^6"^^) = (4.27) 

(Jjdj 

which has the solution 

e{x) = - ^ {3A + 12Ag + 3B{1 + 4^)6"^^ + Ce-^^(A + Be-'^^'f^^) 



^l/4g-3x/2 



3 3 7 
4' 4' 4' 



Ae 



Gx 



B 



(4.28) 



This solution substantially simplifies under the restriction g = —3/4, becoming 



e{x) 



+ 6 



(4.29) 



{B + Ae^^y/^ 

If we choose the integration constant C = 0, then 6{x) = 6, and after elementary 
integration gives (j){x) = \/6x + 0o- For this particular solution, we find the GB 
coupling from (4.4) and ( |4.26[ ) as (up to additive constant) 



i^2(0) 

and the kinetic coupling 



4:^/A 



Fi{4>) 



In 



A + ^A[A + Be 



1/2 



3^ {B + Ae^'t') 



1/2 



(4.30) 



(4.31) 
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where we have set 0o = 0. Thus, there is a particularly simple solution = \/6x + (f)o, 
giving rise to Chaplygin gas cosmology with GB and kinetic couplings as simple 
functions of exponentials of the scalar field. Is worth to mention that this solution 
was obtained in absence of potential, which saves us having to put by hand one more 
input in the model. 

3) The generalized Chaplygin gas solution. 

The generalization of the Chaplygin gas has also been considered to describe in a 
unified way the dark matter and dark energy. It's equation of state is given by 

P=-^ (4.32) 

where p and p are the energy density and pressure of the generalized Chaplygin 
gas, A is a positive constant and a is considered to lie in the range < a < 1, 
which guarantees the stability and causality [51]. The value a = 1 gives the original 



Chaplygin gas. The EoS (4.32) has an equivalent field theory representation in a 
generalization of the Born-Infeld theory [52]. In the scalar field representation used 
in [H2], the Born-Infeld Lagrangian density is reproduced for a = 1. Solving the 
continuity equation and appropriately rescaling the constant A and the integration 
constant B (setting = 1), leads to the Hubble function 

H\x) =[A + 5e-^(°+i)^] ^ (4.33) 



replacing this function in (4.6) gives 



df) 

2{B + Ae^(i+")^) — + 3 (3fi + 4Ae^(^+")^) e{x) + 36 {2gB + (1 + 4^)Ae^(^+")^) = 

(4.34) 
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Integrating this equation gives the following kinetic term 



9{x) 



{B + Ae3(i+")^)2(i+^ 

- 8gB^^ {B + yie^(i+")^)~2(i^ 



1 - 



5 + 2a ^e3(i+°)^- 



2(1 + a)' 2(1 + a)' 2(1 + a)' 



B 



(5 + 2a)55+2 
1 



X 



2^1 



2(1 + a) 



1 + 



2(l + al 



2(1 + a)' 



^g3{l+a)x 

B 



this expression drastically simplifies if g = —3/4 

n„-9x/2 

9{x) 



{B + Ae3("+i)^)2(i+^ 



+ 6 



(4.35) 



(4.36) 



Is remarkable that, again the simplest particular solution to Eq. (4.34) with g = —3/4, 



is 6{x) = 6. In this case the GB coupling, after integration of (4.4) is 

1 1 2 + a Ae3(i+")<^/v^ 



1 + a' 1 + a' 1 + a' 



B 



(4.37) 



and the kinetic coupling from (4.1|) is 

FM = 



1 



_3(1+q), 

A + Be ^ ^ 



(4.38) 



where we used the expression for the scalar field with 00 = 0. The parameters A 
and B can be fixed according to the restrictions imposed by observational data on 



the Chaplygin gas models. Note that in all solutions (4.21), (4.26) and (4.33), the 
universe evolves towards de Sitter phase. Is worth to note that the present reconstruc- 
tion was made without appeal to the potential term, which is an advantage provided 
by the GB coupling. A reconstruction of the Chaplygin and Generalized Chaplygin 
gas cosmologies in scalar field model with kinetic couplings to curvature, have been 
presented in [53] 
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5 Discussion 



We have studied late time dark energy solutions in a scalar field model with kinetic 
terms coupled to curvature and a Gauss Bonnet term coupled to scalar field. The 
GB coupling has the advantage that does not make contributions higher than second 
order (in the metric) to the equations of motion, and therefore does not introduce 
ghost terms into the theory. Given a constant or exponential (of the e-folding vari- 
able) kinetic term, we have obtained explicit solutions for the Hubble function, and 



particularly a new viable solution (4.11 ), leading to different equations of state at dif- 
ferent epochs, including decelerated, accelerated and phantom regimes. A remarkable 



aspect of the solution (4.11) is that we can satisfy any given in advance asymptotic 
and current behavior of the EoS. The model can also be reconstructed for known 
cosmological solutions describing early time radiation (matter) dominance with tran- 
sition to accelerated universe with dark energy dominance, and in some cases even 
with one more transition to phantom phase. This means that we can think the energy- 
momentum tensor as being made of two components: one behaving like a presureless 
matter and other like dark energy with negative pressure. Nevertheless, concerning 
the presureless matter component, in order to be consistent with observations, we 
still need to investigate the clustering properties of the scalar field at scales of the 
order of the cosmic structures we observe in the universe. 

In the e-folding variable we reconstructed the model for a class of important solutions, 
namely for the ACDM-type solution, and for the known Chaplygin and generalized 
Chaplygin gas cosmological solutions. An interesting aspect of these solutions is that 
for the specific numbers g = —3/4 and 6 = 6, the expressions for the scalar field 
and the couplings have been considerably simplified. It is remarkable that in all cos- 
mological scenarios considered, the potential term is absent, which saves us having 
to either implement a restriction on the potential or propose an ad hoc potential. 
Nevertheless, the presence of potential would bring us more possible solutions and 
will be considered elsewhere. 

We have shown that it is possible to construct viable solutions to the dark energy 
problem with different asymptotic behaviors, in the context of the scalar field model 
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with kinetic and GB couplings contributing to the total energy density. All obtained 
solutions for the scalar field and couplings Fi and F2, are simple enough to attract 
the attention on this model. Hence, the present model extends the number of possible 
evolutionary scenarios that explain the nature of the dark energy. 
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